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Abstract. We present HARPS spectroscopy of p. Arae (HD 160691) performed for studying the origin of the metallicity ex- 
cess in this planet-hosting stars. The asteroseismologic campaign led to the previously reported discovery of a 14 earth mass 
planetary companion (Santos et al. 2004b l. The present analysis reinforces this interpretation by excluding other possible pro- 
cesses for explaining the observed Doppler variation and leads to the identification of up to 43 p-mode oscillations with /=0-3, 
frequencies in the range 1.3-2.5 mHz and amplitudes in the range 10-40 cms -1 . 
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1. Introduction 

One particularly important fact that is helping astronomers 
to understand the mechanisms of planetary formation has to 
do with the planet-hosting stars themselves. These stars were 
found to be particularly metal-rich, i.e. they have, on the av- 
erage, a metal content higher, by ~0.2 dex, than the one 
found in field stars without detected planetary companions (e.g. 
Gonzalez [l998l Santos et al. l200T1l2003ll2004ai . 

The impact of this fact is enormous regarding the under- 
standing of the processes of planetary formation and evolution. 
If stars with planets were formed out of more-than-average 
metallic material, then we might conclude that the metallic- 
ity is playing a key role in the formation of giant planets. 
However some evidences suggest that planetary material might 



have fallen onto the planet-hosting stars (Israelian et al. 2001 



Murray & Chaboyer 2002 ) and might be able to produce the 
metallicity excess observed. In this case we might conclude that 
the issue of planetary migration is commonly to be engulfed 
by the star. Which of these two scenarios is correct? In other 
words: is the metallicity excess observed in planet-hosting stars 
of primordial origin, or is it due to stellar pollution effects? 

Asteroseismology consists in measuring properties of p- 
mode oscillations and provides a unique tool for "drilling" stel- 
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* Based on observations collected with the HARPS spectrographs 
at the 3.6-m telescope (La Silla Observatory, ESO, Chile: program 
073.D-0578) 



lar interiors and for testing the internal structure and chemical 
composition of bright solar-type stars (e.g. Bouchy & Carrier 
2003; Bedding & Kjeldsen 2003). It also appears as the only 
means to determine the metallicity gradient across the stellar 
interior in order to solve the question of the origin of metal- 



licity excess in planet-hosting stars (Vauclair 2004, Bazot & 
VauclairHoEb. 



We noted earlier (Bouchy et al.,|2004 Bouchy, 2005) that 
with a radial velocity precision bellow 1 ms _1 , the HARPS 
spectrograph (Pepe et al. 2002 , Mayor et al. |2003 • mounted 
at the 3.6-m telescope offered promising perspectives for such 
asteroseismologic studies. We present here results of our as- 
teroseismologic campaign obtained with this instrument on 
HD 160691, one of the most metal-rich stars with planetary 
companions. Comparison with stellar structure models and in- 
terpretation are presented in a companion paper (Bazot et al. 
2005). 



2. Observations and data reduction 

HD160691 (ju Ara, HR6585) is a G3IV-V star with a V magni- 
tude of 5.1 and a metallicity [Fe/H] of +0.32. A giant planet 
(1.7 Mj up ) in a long period (637 days) orbit was detected 
by Butler et al. (2001). The orbital solution was updated by 
Jones et al. (J2002 1 and Gozdziewski et al. ( 2003 ), who found 
a long term drift possibly due to the presence of a second 
body (> 1.5 Mj up , ~1500 days) in the system. Bright enough 
for HARPS, this target is furthermore well sited for an as- 
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Fig. 1. Radial velocity measurements of HD160691. White cir- 
cles correspond to the night average. The dotted curve corre- 
spond to the orbital fit of the low-mass-planetary companion 
confirmed thanks to additional measurements by Santos et al. 
(2004). 



teroseismology campaign (observable more than 10 hours per 
night). HD160691 was observed over 8 nights in June 2004. 
We took sequences of 100-s exposures with a dead time of 31- 
s in-between. In total 2104 spectra were collected with typical 
night-averaged signal-to-noise ratio (S/N) per pixel in the range 
95-193 at 550 nm. 

The spectra obtained were extracted on-line and in real- 
time using the HARPS pipeline. Wavelength calibration was 
performed with ThAr spectra. The radial velocities were ob- 
tained by weighted cross-correlation with a numerical mask 
constructed from the Sun spectrum atlas. We also computed 
off-line the stellar radial velocities using the optimum weight 
procedure as described by Bouchy et al. (|2001 ) but found no 
significant gain. The journal of the observations is given in 
Table 1 and the radial velocities are presented in Figures^and 
[2] The dispersion of each individual nights, in the range 1.5-2.5 
m s _I , is strongly dominated by the acoustic modes with period 
around 8 minutes as shown in Fig. [3] 

The simultaneous ThAr spectrum was used in order to 
compute the instrumental drift using the optimum weight pro- 
cedure. The spectrograph drift was systematically lower than 
0.75 ms~' during a night and lower than 2 ms _I during the 
whole campaign, illustrating the extreme stability of HARPS. 
Considering this point, we observed without simultaneous 
Thorium during the fourth night in order to construct as more 
properly as possible (without any risk of ThAr contamination) 
a high S/N combined spectrum for spectroscopic analysis of 
HD 160691. 



3. A low-mass-planetary companion 

Our first surprise was to detect in the Doppler signal an un- 
expected modulation with an semi-amplitude of about 4 ms _1 
and a period equal or greater than 8 days. Fig.[2clearly shows 
that we would not have been able to detect such a small signal if 
only one measurement was made each night. Indeed, the semi- 
amplitude of 4 m s corresponds typically to the peak-to-peak 
variation induced by the p-modes in the time series. 
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Fig. 2. Radial velocity time series of individual nights dis- 
played from bottom to top in a relative scale. The box on the 
fourth night is zoomed in Fig. [3] 
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Fig. 3. Zoom of radial velocity measurements showing in the 
time series the presence of p-modes with periods around 8 min- 
utes. The semi-amplitude of about 2 ms~' does not represent 
the individual amplitude of p-modes but comes from the inter- 
ference of several modes. 



In order to confirm the long period modulation, further 
observations were carried on during two GTO runs (program 
073.C-0005) in July and August 2004. Considering that the 
longest periods of p-modes in this star are around 1 1 minutes 
(see Sect.@}, we decided to construct each measurement by av- 
eraging consecutive independent observations over a period of 
twice the longest acoustic period (~ 22 minutes). These further 
measurements allowed us to confirm the modulation and de- 
termine the period to 9.55 days. This Doppler modulation was 
interpreted as caused by a 14 earth mass planetary companion, 
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Table 1. Journal of radial velocity measurements. The mini- 
mum and maximum S/N of spectra obtained during the night is 
given. The dispersion cr RMS is strongly dominated by the seis- 
mic signal and do not reflect the precision of individual point. 



Date 


Nb spectra 


Nb hours 


S/N 


o-rms [ms _1 ] 


2004/06/03 


279 


10.32 


105-181 


2.49 


2004/06/04 


274 


10.19 


65-237 


1.71 


2004/06/05 


285 


10.34 


77-239 


1.79 


2004/06/06 


286 


10.38 


132-235 


2.03 


2004/06/07 


275 


10.35 


43-263 


2.20 


2004/06/08 


229 


8.43 


24-233 


2.04 


2004/06/09 


202 


7.43 


57-147 


1.65 


2004/06/10 


274 


10.42 


32-186 


1.49 



the third one detected in this exoplanetary system and the light- 
est planetary companion ever observed (Santos et al. 2004b I. 

In order to verify if this long period modulation was due to a 
blended binary system or a stellar spot, we conducted a bisector 
analysis of the Cross-Correlation Function (CCF) as described 
in Queloz et al. (2001 1 and Santos et al. (2002i. Such an analy- 
sis is very efficient to discriminate between radial-velocity vari- 
ations due to changes in the spectral-line profiles from pure 
Doppler variations due to the orbital motion of the star. Fig. 0] 
presents a typical CCF of HD 160691. We computed the bi- 
sector velocity of each CCF for 100 different levels (dividing 
the CCF in 100 slices). We defined the "bottom" velocity Vb 
and "top" velocity V, by averaging the values in the range 10- 
35% and 70-95% of the CCF depth, respectively. The differ- 
ence V, - Vb is equivalent to the bisector velocity span and is 
used to measure the variations of the line profile. The result of 
our bisector analysis is presented in Fig. [5] at the same scale 
as Fig [2 The dispersion of the velocity span on an individual 
night is about 2 m s" 1 . No variation appears in phase with the 
radial-velocity variations of period of 9.55 days. We checked 
that no high frequency signals are present in the velocity span 
showing that p-modes induce a pure Doppler variation on the 
spectral lines. The dispersion of the night-averaged-velocity 
span is 0.8 m s _1 . This result indicates that the periodic -radial- 
velocity change is not due to an blended background binary 
(as HD41004 - see Santos et al. 2002 1 nor a dark photospheric 
stellar spot (as HD166435 - see Queloz et al. l2001i . 

We also checked the magnetic activity of HD 160691 by 
computing the chromospheric activity index S based on the 
relative flux level on Call H and K lines (see Fig.|6j. We mea- 
sured the flux in two 1 A pass bands centered on the H and K 
lines normalized by two 20 A wide sections of photospheric 
flux centered at 3900 and 4000 A. The index S , computed on 
the 2104 available spectra, has values in the range 0.120-0.125. 
Once in the Mount Wilson scale (Vaughan et al. 1978 1, this cor- 
responds to a logfl^=-5.034±0.006 (Noyes et al. [1984! ). As 
mentioned by Santos et al. (2004b I this value, in close agree- 
ment with the one obtained by Henry et al. ( 1996 1, is typical of 
an inactive solar type star. 
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Fig. 4. Typical Cross Correlation Function of HD160691 with 
location of the intervals used to compute the bisector velocity 
span V t - Vb- 




Fig. 5. Bisector velocity span of the CCF profile of HD 1 6069 1 . 
White circles correspond to the night average. 




3964 3966 3968 3970 3972 




, 1 , , , 1 , , , 1 , , , 1 , , , 1 U 

3930 3932 3934 3936 3938 

Wavelength [A] 

Fig. 6. Call H and K lines for HD 160691 showing no chromo- 
spheric activity. 
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Fig. 7. Power spectrum of radial velocity measurements of 
HD160691. 

Finally our acoustic spectrum analysis (see Sect. |4ji allows 
to estimate the rotational period at ~22 days which definitively 
eliminates any stellar activity origin in the 9.55 days signal. 

We noticed that the residuals around the best Keplerian fit 
have a rms of 0.43 ms 4 for the 8 night-averaged measure- 
ments made during our campaign of asteroseismology (aver- 
age of more than 200 individual observations) and have a rms 
of 1.3 ms~' for the other nights (average of 15 observations). 
This is due to the presence of low frequency (few hours) modu- 
lations in the Doppler signal (see Fig.|5J with semi-amplitudes 
of 1-2 ms~' which are not averaged by 15 consecutive inde- 
pendent observations. We checked that this noise did not ap- 
pear in the Thorium signal. We also checked that this noise was 
not introduced by chromatic effect on the fiber entrance due 
to atmospheric-dispersion-correction residuals with change of 
seeing or change of fiber centering. The flux ratio between blue 
and red spectral orders do not present correlation with the low 
frequency modulation seen in the Doppler signal. The origin of 
this noise is still unclear, but considering that it does not seem 
to be an instrumental effect and that its amplitude change with 
time, we strongly suspect that it could have a stellar origin (like 
granulation noise). 

A third exoplanet with low-mass (14 M B ) and short-period 
(9.55 days) was not completely unexpected and was considered 
as a possible lucky by-product of this asteroseismology cam- 
paign. In fact statistics of the exoplanet sample show that the 
probability to find a new exoplanet is higher around a planet- 
hosting star (more than 10% of known exoplanets are in a mul- 
tiple system), even if the host star is overmetallic. We show 
here that the only way to find such an companion with a semi- 
amplitude of only 4.1 ms _1 was to perform such an intensive 
campaign of radial velocity measurements in order to average 
properly and completely the stellar oscillations. 

4. Acoustic spectrum analysis 

In order to compute the power spectrum of the velocity time 
series of Fig. ^ (corrected from the orbital motion), we used 
the Lomb-Scargle modified algorithm (Lomb 119761 Scargle 
I1982 i for unevenly spaced data. The resulting LS periodogram, 
shown in Fig.0 exhibits a series of peaks between 1.3 and 2.5 
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Fig. 8. Sum of the echelle diagram showing the modes 
/=0, 1,2,3 and their side-lobes. 

mHz modulated by a broad envelope, which is the typical sig- 
nature of solar-like oscillations. This signature also appears in 
the power spectrum of each individual nights. The mean white 
noise level computed in the range 0.7-1.2 mHz is 2.6 x 10~ 3 
m 2 s -2 . Considering a Gaussian noise, the mean noise level 
in the amplitude spectrum is 4.5 cms -1 . Considering that the 
time series is based on 2104 measurements, the corresponding 
velocity accuracy thus corresponds to 1.17 ms _1 . The uncer- 
tainty due to photon noise is estimated to 0.52 m s indicating 
that we are not photon-noise limited. The origin of the external 
noise level of ~1 ms _1 is still unclear. The main instrumental 
limitation of HARPS could comes from the guiding and cen- 
tering errors on the optical fiber input but tests made during 
commissioning indicated error sources at the level of 0.2 m s" 1 
(Mayor et al., 2003). We strongly suspect a stellar origin of this 
noise (like granulation noise). 

In solar-like stars, p-mode oscillations are expected to pro- 
duce a characteristic comb-like structure in the power spectrum 
with mode frequencies v„j,„, reasonably well approximated by 
the simplified asymptotic relation (Tassoul 1980): 

/ 1(1+1) 

v n j, m * Av • (n + - + e) - 6v 2 7 ^ m (!) 

2 6 

The two quantum numbers n and / correspond to the radial 
order and the angular degree of the modes, respectively. For 
stars of which the disk is not resolved, only the lowest-degree 
modes (/ < 3) can be detected. In case of stellar rotation, p- 
modes need to be characterized by a third quantum number m 
called the azimuthal order (— / < m < I). Quantities Av, 6vq2 
and Q correspond to the large spacing, the small spacing and 
the rotational splitting, respectively. 

The power spectrum, shown in Fig.0 presents a clear and 
unambiguous periodicity of 90 fiHz both on its autocorrelation 
or in the comb response. In order to identify the angular de- 
gree I of each mode individually, we cut the power spectrum 
into slices of 90 //Hz and summed up them. The result, which 
corresponds to the sum of the echelle diagram, is presented on 
Fig.|8]and allows to identify unambiguously modes / = 0, 1, 2, 3 
and their side-lobes due to the daily aliases at +1 1.57 yiiHz. 

We selected the strongest peaks with a signal-to-noise ratio 
greater than 3 in the amplitude spectrum and identified them 
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Table 2. P-mode frequencies for HD160691 (jxHz). The fre- 
quency resolution of the time series is 1 .56 //Hz. Frequencies 
with label a correspond to peaks with a signal-to-noise ratio in 
the range 2.5-3 and have to be considered with caution. 



Table 3. P-mode amplitudes for HD160691 (cms 1 ). 
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Fig. 9. Power spectrum of HD 160691 with the identified fre- 
quencies in Table [2]marked by dotted lines. 
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as p-modes with n value based on the asymptotic relation (see 
Eq. (1)) assuming that the parameter e is near the solar value 
(e ~ 1.5). The Frequency and amplitude of these 33 p-modes 
are listed in Tables|2]and[3] The amplitude was determined by 
assuming none of the p-modes are resolved and corresponds to 
the height of the peak in the power spectrum after quadratic 
subtraction of the mean noise level. Considering that the fre- 
quency resolution of our time series is 1.56 //Hz, we adopted an 
uncertainty on no-resolved oscillation modes of 0.78 //Hz. Such 
an uncertainty can be considered as conservative in case of in- 
finite time-life modes with high signal-to-noise ratio (>3.5). 
A second inspection of the power spectrum with selection of 
peaks with signal-to-noise ratio in the range 2.5-3 (with am- 
plitude in the range 11-13.5 cms 4 ) allowed us to propose 10 
additional p-modes. We are aware that some of these weaker 
modes could be miss-identifications and should be considered 
with caution and lower confidence. 



10 20 30 40 50 60 70 80 90 

Frequency modulo 90 /j-Hz [^iHz] 

Fig. 10. Identified p-modes 1=0 (□), 1=1 (■), 1=2 (a), and 1=3 
(a) in echelle diagram format. The size of the symbols are pro- 
portional to the amplitude of the modes. Crosses correspond to 
peaks not identified in the power spectrum. 

All the identified modes are displayed on the power spec- 
trum in Fig. [5] and in the echelle diagram in Fig. ^3 A few 
peaks with signal-to-noise greater than 3 but not identified as 
p-modes or side-lobes are displayed in Fig. EH as crosses. It 
could indicate that some modes are split due to their damping 
time. 

The large spacings Av = v„+\j - v n j were computed for 
each pair of consecutive modes with same angular degree. For 
the 1=2 modes, which are split into 2 components due to the 
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Fig. 11. Large spacing Av versus frequency for p-modes of de- 
gree Z=0 (□), 1= 1 (■), 1=2 (A) and 1=3 (A). 

stellar rotation, we used the averaged value of the 2 frequen- 
cies. The result is displayed in Fig. ^2 Neglecting the first three 
points, the residual of the large separations after a third order 
polynomial fit is 1.04 //Hz RMS. This value indicates that the 
adopted uncertainty of 0.78 x V2 = 1.1 /iHz for the large sep- 
aration is appropriate. The small spacings 6vo,i = v„+i,o - v n ,2 
and (5vi,3 = v n +i,i - v n ,3 are displayed in Fig.[^] The average 
small spacings 8vq2 and 8v\$ are 5.7 + 0.4 //Hz and 7.5 ± 0.4 
//Hz respectively. 

The splitting of 1=2 modes are displayed in Fig. [Q] 
Considering that such a splitting does not appear in modes 
1=1 and the fact that the visibilities of modes to = +2 are 
higher than modes to + 1, we strongly suspect that the observed 
1=2 modes correspond to azimuthal order m = +2. The av- 
erage splitting of 2.1 ± 0.2 //Hz leads to a rotational splitting 
Q. = 0.53 + 0.05 //Hz which corresponds to a rotational period 
of 22 + 2 days in agreement with the very low magnetic activity 
found in Sect. Modes /=3 do not appear to be split. It could 
indicate that we only observe modes (I = 3, to = 0). An alter- 
native is that we observe modes (I = 3, to = -3). In this case 
modes (/ = 3, to = +3) are expected to be at — 11.1 //Hz from 
the / = 1 modes, hence completely merged with the /= 1 side- 
lobes. 

5. Conclusions 

Our spectroscopic observations of HD 160691 first yield to the 
discovery of a very light planetary companion (14 M @ ) with 
short period (9.55 days) (Santos et al. 2004b I. Our analysis en- 
forces the validity of this planetary companion excluding all 
sources of ambiguities. Our observations of HD160691 also 
yield a clear detection of p-mode oscillations. Up to 43 p- 
modes have been identified in the power spectrum between 1.3 
and 2.5 mHz with an average large spacing of 90 //Hz, an aver- 
age small spacing of 5.7 //Hz, a rotational splitting of 0.53 //Hz 
and an envelope amplitude of about 30 cms 1 . The identified 
p-mode frequencies will be compared with stellar models in a 
companion paper (Bazot et al. 12005b . 

Acknowledgements. The authors wish to thank R.L. Gilliland for his 
helpful suggestions. We are grateful to ESO staff support at the 3.6- 
m telescope. S.V. acknowledges a grant from Institut Universitaire 




— 4 - 

i 

' 3 - 
tafl 

I 2 
ft 1 - 

in 

o — 

1.2 



1.6 1.8 2 2.2 

Frequency [mllz] 



Fig. 12. (Top panel) Small spacings (5vn,2 (•) and 6v\ 3 (o) ver- 
sus frequencies. (Bottom panel) Splitting of modes 1=2, m=+2 

(A). 



de France. N.C.S. would like to thank the support from the Swiss 
National Science Foundation and the Portuguese Fundacao para a 
Ciencia e Tecnologia in the form of a scholarship. 



References 

Bazot, M., & Vauclair, S„ 2004, A&A, 427, 965 

Bazot, M., Vauclair, S., Bouchy, F, et al., 2005, A&A, in press 

Bedding, T, & Kjeldsen, H., 2003 , PASA, 20, 203 

Bouchy, F, Pepe, F, Queloz, D., 2001, A&A, 374, 733 

Bouchy, F, & Carrier, F, 2002, A&A, 390, 205 

Bouchy, F, & Carrier, F, 2003, Ap&SS, 284, 21 

Bouchy, F, Maeder, A., Mayor, M., et al, 2004, Nature, 432(7015), 2 

Bouchy, F, 2005, in Dome C Astronomy/Astrophysics Meeting, EDP 

Sciences, EAS Publications Series, Vol.14, in press 
Butler, R.P., Tinney, CO., Marcy, G.W., et al., 2001, ApJ, 555, 410 
Gonzalez, G., 1998, A&A, 334, 221 

Gozdziewski, K., Konacki, M., & Maciejewski, A.J., 2003, ApJ, 594, 
1019 

Henry, T.J., Soderblom, D.R., Donahue, R.A., et al., 1996, AJ, 111, 
439 

Israelian, G., Santos, N.C., Mayor, M., et al., 2001, Nature, 411, 163 
Jones, H.R.A., Butler, R.P., Marcy, G.W., et al., 2002, MNRAS, 337, 
1170 

Lomb, N. R., 1976, Ap&SS, 39, 447 

Mayor, M., Pepe, F, Queloz, D., et al., 2003, Messenger, 114, 20 
Murray, N., & Chaboyer, B., 2002, ApJ, 566, 442 
Noyes, R.W., Hartmann, L.W., Baliunas, S.L., et al., 1984, ApJ, 279, 
763 

Pepe, F, Mayor, M., Rupprecht, G., 2002, Messenger, 110, 9 
Queloz, D., Henry, G.W., Sivan, J.P., et al., 2001, A&A, 379, 279 
Santos, N.C., Israelian, G., & Mayor, M., 2001, A&A, 373, 1019 
Santos, N.C., Mayor, M., Naef, D., et al., 2002, A&A, 392, 215 
Santos, N.C., Israelian, G., Mayor, M., et al. 2003, A&A, 398, 363 
Santos, N.C., Israelian, G., Mayor, M., et al. 2004a, A&A, 415, 1153 
Santos, N.C., Bouchy, F, Mayor, M., et al., 2004b, A&A, 426, L19 



Bouchy et al.: The acoustic spectrum of /i Arae 



7 



Scargle, J. D., 1982, ApJ, 263, 835 
Tassoul, M, 1980, ApJS, 43, 469 
Vauclair, S., 2004, ApJ, 605, 874 

Vaughan, A.H., Preston, G.W., & Wilson, O.C., 1978, PASP, 90, 267 



